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We derive analytic expressions for the possibility, probability, and expected size of global spread- 
ing events starting from a single infected seed for a broad collection of contagion processes acting 
on random networks with both directed and undirected edges and arbitrary degree-degree correla- 
tions. Our work extends previous theoretical developments for the undirected case, and we provide 
numerical support for our findings by investigating an example class of networks for which we are 
able to obtain closed-form expressions. 
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I. INTRODUCTION 

Spreading mechanisms playing out on generalized ran- 
dom networks constitute a rich and compelling class of 
tractable contagion models [UH]. First, while real world 
complex networks are rarely, if ever, pure Erdos-Renyi 
networks, they often possess a strong, describable mea- 
sure of randomness 3J, once the dominant aspect of 
degree distribution is acknowledged Second, simple 
models of network-based spreading have yielded impor- 
tant insights into spreading phenomena such as the 
spread of infectious diseases [5] [5] , cascading failures in 
power grids [5], and social contagion processes 0- 
118] . Finally, many random network models are amenable 
to analytic investigations and researchers have natural- 
ly built on areas of statistical mechanics — with its great 
tradition of exactly solvable models — such as the study 
of percolation on lattices |19j . 

Here, we examine contagion processes acting on mixed 
directed and undirected degree-assortative random net- 
works. Specifically, for the case of a single seed, we derive 
and verify by simulations analytic expressions for three 
key aspects of these systems: (1) the possibility of a glob- 
al spreading event; (2) the probability of a global spread- 
ing event; and (3) the expected final size of a successful 
global spreading event. We make the distinction between 
possibility and probability, the former referring to the 
potential for spreading (i.e., whether or not the system is 
in a phase where spreading may occur), and the latter to 
the quantified chance that a macroscopic spreading event 
may arise given the nature of the initial seed (e.g., ran- 
dom or targeted). To put it another way, possibility is a 
categorical yes/no criterion and probability is a quantita- 
tive one; they ask different kinds of questions, and elicit 
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different kinds of analyses for their determination. Thus, 
while we could simply derive the probability of global 
spreading only and thereby immediately know if global 
spreading was possible or not (corresponding to non-zero 
or zero probabilities), obtaining the possibility of global 
spreading alone is important as it directly reveals phase 
transitions, and further involves a transparent, physically 
argued calculation 20 . 



We base our work most strongly on two groups of 
authors' findings: Boguha and Serrano [21 , who pro- 
vided a general formulation for such networks; and Glee- 
son and Cahalane [23J , who derived the final size of 
global spreading events for general contagion models on a 
wide array of network structures, including the social-like 
threshold model on random networks [T2J \T5\ [23] . Out- 
work is related to that of Meyers et al. [25] who examined 
disease-spreading models on mixed directed and undi- 
rected uncorrelated networks; our analytic methods are 
essentially disjoint and we treat more general spreading 
mechanisms (e.g., social- like ones), while Meyers et al. 
explored various real-world applications. 



We structure our paper as follows. In Sec. [IT] we define 
the family of random networks and contagion processes 
we investigate here; in Sec. |III| we provide physically- 
motivated expressions for the possibility and probability 
of global spreading events starting from a single seed; 
and in Sec. |IV| we derive coupled evolution equations 
that describe the growth of a global spreading event, as 
well as yield the expected final size. In the appropriate 
limits, our equations collapse to those for various net- 
work subclasses involving purely directed or undirected 
links, and correlated or uncorrelated nodes. In Sec. [V] 
we obtain exact results regarding spreading on a specific 
family of random networks. We close with a few remarks 
in Sec. Ell 
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II. MODEL DESCRIPTION 
A. Generalized random networks 

We consider random networks containing undirect- 
ed and directed links along with arbitrary correlations 
between nodes based on degrees. Following Boguna and 
Serrano [21j . we allow each node to have k u undirect- 
ed edges, fcj incoming directed edges, and k Q outgoing 
directed edges. We assume all edges are unweighted. 
We represent a node by the generalized degree vector 
k = [k u , ki, k ] T , which we will refer to simply as a node's 
degree, and we write P(k) for the degree distribution. 

To account for correlations, three conditional probabil- 
ities are needed: P^(k | k'), P^(k\k'), and P^(k\k'); 
these quantities give the chances that an edge starting at 
a degree k' node ends at a degree k node and is respec- 
tively undirected, incoming, or outgoing relative to the 
destination degree k node (note that this convention for 
directed edges is opposite that used in [2"T]). 

For networks to be well defined (i.e., realizable), these 
probabilities must be constrained by two detailed bal- 
ance equations. In determining the probability that an 
edge of a certain type runs between nodes of degree k 
and k' , we must obtain the same result whether we start 
at the former or latter node. We write the probability 
that a randomly chosen edge is undirected and connects 
a degree k and degree k' node as P( n '(fc, k'). Noting that 
the probability that a random end of a randomly selected 
undirected edge emanates from a degree k node is given 
by kn ifc P ; we first have that 

P W(k,k>) = pW(jf|j?)^p 

= p(»>(#|jfe)*^^ = p<»)(je>,jfe). (i) 

(feu) 

For directed edges, we define p( dir )(fc,fc') as the proba- 
bility that a randomly chosen edge is directed and leads 
from a degree k' node to a degree k node. Similar to 
the balance equation for undirected edges, we use the 



quantities 



fc p(fc) 



and 



fciP(fc) 



which give the probabilities 



<fe ) — <fe,; 

that in starting at a random end of a randomly selected 
edge, we begin at a degree k node and then find ourselves 
travelling (1) along an outgoing edge or (2) against the 
direction of an incoming edge. We therefore have 



p( dir )(fc,fc') = P«(fc|fc') 



H k' P(k') 



p( o) (k'\k) 



k { P{k) 



(K) v 1 ' (h) ' 

(2) 

Note that since (fc u ) = (k' u ) and (k' Q ) = (k Q ) = (ki), 
the denominators in Eqs. ([T]) and ([2| are equal and may 
be omitted [3T]. Further, our alternate definitions of 
pW(fc' | jfe) and P(°)(fc' | k) mean that Eq. {i| has a form 
different to that given in [2"T] . 



For the class of random networks given above, Boguna 
and Serrano determine a number of structural results 
regarding percolation, including the sizes of the giant 
in-component, out-component, and strongly connected 
component |21j . Our goal here is to examine the behav- 
ior of generalized spreading processes on such networks, 
and we describe these next. 



B. Contagion processes 

We consider synchronous discrete time contagion pro- 
cesses, though our results can at least in part be extended 
to asynchronous discrete and continuous time process- 
es [23 . We assume that once nodes are infected, they 
remain so permanently, an aspect that is needed for com- 
puting the final size of a global spreading event. We write 
the probability of node j becoming infected in time step 
t + 1 as 



Bj\ki n {] k u + ki), 



(3) 



given that fc; n f of node j's total of k u + ki undirect- 
ed and incoming edges emanate from infected nodes at 
time t. Here, Bj is an arbitrary, node-specific 'response 
function' mapping to the unit interval. Now, for the 
general class of contagion models we consider here on 
infinite random networks, we need to know only the 
average response function for each node subclass. Tak- 
ing all nodes of degree k, having indices in the set 

J k = {jfc,i'Jfc,2' • ■ ■ >h,n' ■ ■ ■}' we com P ute tnis average 
response function as 



B h nt ,% = - BjihnCK + h) 



(4) 



The quantity B k f g is then the probability that a ran- 
domly chosen node of degree k is infected at time t + 1 
given that at time t, it has fci n f infected incoming and 
undirected edges. 



III. POSSIBILITY AND PROBABILITY OF 
GLOBAL SPREADING 

In (20] . we derived a global spreading condition for 
discrete and continuous time contagion processes with 
the possibility of recovery acting on generalized random 
networks. Defining a = (y, A) to represent a pairing 
of a type v node and type A edge, we argued that the 
number of infected node-edge pairs fs grows as a func- 
tion of network distance d from a seed as fs(d + 1) = 
^2s' P-aS' fa' (d) , where R&s> depends simply on network 
structure and the spreading process [20]. As a special 
but still broad case, we showed that for the networks we 
consider here, the growth rate equation for the number of 
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infected edges emanating from degree k nodes a distance 
d from an initiating node obeys the following: 



4 u) ( rf - 



- R Sfe' 



4 n) (d) 



(5) 



where 



(6) 



p( u )(fc|fc').(A: u _ 
p( u )(fc|fc')»fco 



P«(fc|fc') 



Here, the quantities /^(d) and f^\d) are the num- 
ber of 'infected' undirected and outgoing edges leaving 
an infected degree k node a distance d steps from the 
seed. We have expressed the form of Rj£, so as to make 
clear the three components making up general spreading 
conditions: (1) probability of connection [p( u )(fc | k') and 
PW(fc | k')]; (2) resultant newly infected edges [(fc u — 1), 
k u and k a factors]; and (3) the probability of infection 
{B 1 g) [10]. The above agrees with the contagion condi- 
tion found earlier by Boguha and Serrano for the emer- 
gence of the giant out-component using a generating 
function approach. Note that these calculations depend 
on the local pure branching structure of random networks 
with zero clustering; for recent advances for the non-zero 
clustering case, see [2"rJH2"9"] . 

The full gain matrix R and edge infection counts f(d) 
can be laid out as follows: 



R 



^Tsifei '^ J kik2 



and f(d) 



f { i\d) 

fhd) 



(7) 

The condition for the possibility of global spread- 
ing events is therefore that the maximum eigenvalue of 
[Rpt/1 exceeds 1: 



sup{H : I* € <r ([*&])} > 1 



(8) 



where cr(-) indicates spectrum. 

Next, we determine the probability of a global spread- 
ing event given the initial seed is of degree fc and hence 
the overall probability given a randomly selected seed; 
we refer to these quantities as 'triggering' probabilities. 
While in determining the probability of a global spread- 
ing event we must also determine the possibility, the 
direct calculation we have just presented for the latter 
is needed to demonstrate a physically-motivated clarity. 

We define to be the probability that an infected 

undirected edge leaving a degree fc node will lead to a 
giant component of infected nodes. Similarly, we define 



to be the probability that an infected outgoing edge 

from a degree fc node will generate a global spreading 
event. Using the Markov nature of random networks, we 
can write down recursive, closed-form relationships for 
these two probabilities: 



(u) 



E p(u) (fc'ifc) 



)(°) 



and 



Qf 



J2 pii) (k'\k) 



1 - 



(9) 
B il» 



(10) 



In these equations, we have encoded the understand- 
ing that if an infected edge generates a global spread- 
ing event, then it must infect its target node which in 
turn m ust be successful in infecting its other neighbors. 
In Eq. ( 10 1 , for example, pW (k'\k) is the probability that 
the undirected edge leads from an infected degree k node 
to a degree k' node which it infects with probability B 1 p . 

The quantity (1 - Q { p) k '"(l - Q ( p) k '° is the probability 
that none of the infected node's other undirected or out- 
going edges successfully spread the infection, and hence 



( u )\fc' 



(1 



is the probability that at 



least one does. 



Both Q ( - u) 

k 



and Q\ 

k 



(o) 



can be determined from Eqs. (9) 
and ( 10 1 cither numerically or exactly (as per our exam- 
ple later in Scc|v]). Having done so, we can then compute 
the probability that infecting a single degree k node trig- 
gers a global spreading event: 



Qf) k ° 



(11) 



which is the complement of (1 — Q^) ku (l 



4 0) ) fe °' 

the probability of failure to trigger. The probability 
that infecting a randomly chosen node triggers a glob- 
al spreading event is then simply Prig = 2g Ptrlg(k), 
or 



p tris = J2 p (k) 



(12) 



In similar fashion, the triggering probability for non- 
random, strategic selections of the initial seed can read- 
ily be obtained. Appropriate limits of Eq. ( 12 ) also 



recover triggering probabilities for simpler families of ran- 
dom networks such as undirected, uncorrelated networks 
with prescribed degree distributions. Finally, considering 
the limit of Ptrig - ► retrieves the condition for global 
spreading found above. 
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IV. FINAL SIZE OF SUCCESSFUL GLOBAL 
SPREADING EVENTS 



as 



We complete our main analysis by determining the 
final size of a global spreading event building on the work 
of Gleeson and Cahalane [55] and later Gleeson [53]. We 
shift our focus from spreading away from a seed (expan- 
sion) to spreading reaching a node (contraction). 

We consider an arbitrary fixed node in the network 
and compute the probability that incoming edges (direct- 
ed or undirected) are infected and sufficient in number 
that the node itself becomes infected at a certain time. 
To do so, we need to first determine the probabilities 
that undirected and incoming edges arriving at a degree 
k node are infected at time t, 9^ and 9~ . As with the 

' k,t k,t 

possibility and probability of spreading, edge-edge tran- 
sitions are the best framing for this calculation. Edges 
will be infected at time t + 1 if the node from which they 
emanate becomes infected in that time step, and this 
in turn depends on the infection levels of the incoming 
edges. Assuming a fraction 4> > of initially infected 
seeds in the network, we obtain the following expression 
for the fraction of infected directed and incoming edges 
in the network at time t + 1: 



a EE (" u )(": 



k,t 



Ju 



1 



fc,tj L mJ I 



(fe-ii) 
x Bju+juK+h- (15) 



9 (i) 
k,t 



q(u) _ 



To determine the final size, we set 9\ 

? « 
fe',t+l 



= 

k't 



)(u) 
/c' ,oo 



and 



9 (0 

fc'.c 



solve for the steady-state solutions #l u ^ 

fc' ,oo 



in Eqs. (|13|) and (|14|) and 
Sub- 



and # 

k' ,oo 



stituting these values into Eq. ( 15 ) gives us the expected 
final size </>oo which is, among other things, a function of 
cf>o, the fraction of nodes initially infected. For the single 
seed case we consider in this present work, the final step 
therefore is to take the limit <fio —> 0. Note that as for the 
triggering probability, the condition for global spreading, 
Eq. ( Ts|) can be recovered by linearizing Eqs. (13), (14 1, 
and |l5]) (see ref. [53J). 



9 (u) 
fc,t+i 



(i-^E^^i^EEf^jf. 





'1 


1-^1 






fc'.tj 



(fci'-iO 



and 



9 « 
fe,t+i 



d-«E^(^)EE(. : )(- : 























fc/.t. 





x B Ju+j . ifci+i . ; . (14) 

Since we are now considering contraction rather than 
expansion, more than one edge may contribute to the 
infection of a node, hence the sum over nearly the full 
range of infection probabilities, the {B^+^^+k'.}- 

The overall fraction of infected nodes at time t, equiv- 
alently the probability that a randomly chosen node 
becomes infected at time t, depends on and 9y 

' 1 k't k't 



V. EXACT SOLUTION FOR AN EXAMPLE 
DEGREE-CORRELATED RANDOM NETWORK 
WITH MIXED DIRECTED AND UNDIRECTED 
EDGES 



To test our analytic expressions for the possibility, 
probability, and expected final size of a global spreading 
event, we consider a family of general random networks 
for which our equations are exactly solvable. As shown 
schematically in the margins of Fig. [T] we allow four types 
of nodes with the following degree vectors, which, again, 
have the form k = [k u , ki, k ] T : 



"2" 




"o" 




"o" 


1 


, k 2 = 







1 


1 




1 








ki 



and which occur with abundances 



and kii = 



(16) 



P(ki) = \,P(h) = l,P(h) - \, and P(k 4 ) = | 



(17) 



We define the degree-degree conditional probabilities as 
dependent on two tunable parameters, r un( j and r^ r : 



P {u) (k\k') 



Tund 0(1- T und )' 





.(1 - Tund) T und . 



, (18) 
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FIG. 1: Central plot: For the toy network model described in Sec. [V] final size 4>oo as a function of the model parameters 
T un d and Tdir- Size is mapped to a linear gray scale with white indicating no global spreading. The dashed line marks the 
theoretically determined phase transition given in Eq. ( |22[ ). The example networks shown around the plot give a sense of 
the kinds of networks realized for the corresponding r un d and Tdir- Note that we show exact forms for networks with 20 
nodes to make clear how node types are correlated in phase space, and forms will differ for larger networks. For example, for 
T" U nd = Tdi r = 1, networks will comprise a giant component of type 1 nodes (open circles) with the remaining three node types 
represented only in isolated pairs. Simulation details: We formed each network with N = 10 4 nodes made up of a 1:1:1:2 ratio 
of node types 1 through 4. We constructed each network initially to have r un d = Tdi r = 1, which was simple algorithmically, 
and then shuffled edges until desired values of r U nd and Tdir were reached, using an approach similar to those described in |30| 
and [3T]. We further shuffled each edge type 10,000 times to ensure randomization. For each r un d and Tdi r in 0, 0.01, 0.02, . . . , 
1.00, we generated 100 networks and randomly picked 1000 seeds for a total of 10 s samples. 



and 



P^(k\k')] = \p {o \k\k') 



Tdir 


(1 - T dir ) 




















- Tdir) 


Tdir 





















(19) 



where < T^ ir , r un( j < 1, and k and k' correspond to rows 
and columns. 



We have chosen Tdir and r un d so that increasing them 
will tend to increase global connectivity, with r un d con- 
trolling correlations between nodes through undirected 
edges, and through directed ones. There are four 
clear limiting cases, as shown in the corners of Fig. [T] 
For example, when rair = T un d = 1 (upper right corner of 
Fig. [I]), type 1 nodes are connected only to other type 1 
nodes creating a giant component, while the other three 
types combine to form isolated pairs with either directed 
or undirected connections. At the other extreme when 
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Tdir = r und = 0, (lower left corner of Fig. [I]), each of 
the four edges from type 1 nodes connect only to type 
2, 3, and 4 nodes, meaning the network is composed of 
discrete, five-node components. The six other example 
networks in Fig. [I] give a sense of the other possible con- 
figurations contained within this simple network family 
we have constructed. 

We obtain results for general response functions, while 
for comparison with simulations, we consider a test 
contagion process with the following single-parameter 
threshold transmission probabilities: 

B Q g. = 0, B x ^ = ft, and B. %.=! otherwise. (20) 

where i = 1,...,4. The choice B Q g =0 means no 
nodes spontaneously become infected (as might model 
the action of an exogenous source of infection). In the 
case that [3 = 1, then this set of responses means that 
if a node finds at least one neighbor at the end of an 
undirected or incoming edge that is infected, then the 
node itself becomes infected in the next time step. For 
ft < 1, a random fraction ft of degree k\ nodes become 



infected in the time step following the infection of a single 
neighbor whereas 1 — ft remain uninfected. As discussed 
in Sec. |IIB[ individual response functions need only give 
this average response function; for example, a fraction 
ft of degree ki nodes might have a deterministic thresh- 
old of 1 with the remaining fraction of 1 — ft having a 
deterministic threshold of 2. 

Returning to Fig. [T] the gray-scale plot shows the frac- 
tional size of successful global spreading events as a func- 
tion of T un d and Tdir for the specific spreading mechanism 
described above with ft(= B 1 r ) = 1. We see a clear 
phase transition indicated by the dashed curve and our 
next task is to find its analytic form. 



A. Global spreading condition 

Using the spreading conditions contained in Eqs. (|5]), 
([7]), and ([8]), and the model's definition, we find that glob- 
al spreading may occur when the maximum eigenvalue of 
the following gain matrix exceeds unity: 



R- - R- - R- 

fcifci h\h-2 fc'ifc'3 



R 



^£2 ki "^£12 £2 



hi &4 
k 2 £4 
^fe 3 £4 

^^4^4 



2T dit B lii 2(1 - Tdir)^ g 2 {l-T nnd )B 1 ^ 0" 
Tun d B irki r dil B lA (l-r dlr )S a2 (l-r und )B a4 
00 00 







Clearly, only the top left hand corner of this gain matrix 
matters as global spreading, if possible, must occur on 
a giant component. Upon substitution of the model's 
response functions given in Eq. (20), we find the largest 



eigenvalue is \ (r und + r dir + y/ (r und + r dir ) 2 + 4r und T dir )/3 
and we find that global spreading events therefore occur 
in the region described by 



(l + T und ft)(l + T dil ft) > 2. 



(22) 



For the case ft(= B x g ) = 1, this equation is indeed rep- 
resented by the dashed curve shown in Fig. [l] perfectly 
matching the phase transition demonstrated by our sim- 
ulations. We can also now readily determine that spread- 
ing may occur for some values of r und and r d ; r providing 
0>y/2-l. 



B. Probability of global spreading 

In computing the probability that a degree k node ini- 
tiates a global spreading event, we observe that because 
only type 1 nodes can transmit an infection, we need 
only solve the recursion equations given in (|9| and ( 10 ) 



for Qf 1 ' and . Nodes of type 2 and 4, possessing one 

fci k± 

outgoing and one undirected edge respectively, may trig- 
ger global spreading but obviously cannot be involved in 
transmission, and nodes of type 3 can neither start nor 
help spread an outbreak. Eqs. ([9| and ( 10 ) reduce to the 
nonlinear coupled equations: 



} (u) 



Tund 



and 



Qf = r dir 



1-1- 



)(o) 

'fcl 



(°) 



The equations are solvable and we find 



'fci 



1 



7 7un d / 



with 



^ \Tundft 



- 1 



1 , fl x2 T und 1 
[Tundft) + 

4 Tdir T dir p 



(23) 



ft. (24) 



(25) 



(26) 
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For (3 = r und = r dir = 1, 



)(u) 



1. In turn, Q { ° ] 

k2 



and can be expressed in terms of and Q^: 

&4 fei Ail 



Q(°) = (l-r dir ) 



and 



= (l - r und ) 



1 - 



)(o) 



)(o) 



P, (27) 



/?. (28) 



A first check on these triggering probability expressions 
is that they are in agreement with the phase transition 



recorded in Eq. (|22l): in other words, and should 

k i k i 



vanish along the phase transition. We see that upon set- 
ting the right hand side of Eq. (251 to zero, rearrange- 
ment indeed leads to the condition(l+T un d/9)(l+rdir/5) = 
2. 

We compute the triggering probability given a random- 
ly chosen seed using Eq. (12): 



trig 



1 



1 



)(u) 



1 



ki 



5^ 



(o) 



(u) 



(29) 



We compare our theoretical computation of Ptrig with 
simulations in Fig.[2]for one transect in r und — r d i r param- 
eter space (T d i r = 0.66, r und varying), and some example 
values of j3. 

Note that for (3 = 1, all nodes are vulnerable provid- 
ing they can be reached along an edge, and the conta- 
gion model's behavior exhibits an inherent symmetry in 
the network's giant in-component and out-component, 
implying that -Ptrig = <j>oo- For the infection probabili- 
ty Ptrig) the initial node must be part of the giant in- 
component which is made up of type 1, 2, and 4 nodes. 
The hnal infected component will match the giant out- 
component which is turn made up of type 1, 3, and 4 
nodes. The giant strongly connected component is found 
in the intersection: type 1 and 4 nodes. 




0.4 0.6 
7"und 



FIG. 2: For the model described in Sec. [V] the probabil- 
ity that infecting a randomly chosen node leads to a glob- 
al spreading event, Ptrig, as a function of undirected edge 
assortativity r un( j, a fixed value of directed edge assortativi- 
ty, Tdi r = 0.66, and varying values of j3 = B l g . The curves 
correspond to outp ut from simulations [squares] and theory 
[solid line, Eq. |29|]. For the simulation results, we recorded 
a successful global spreading event if the final size exceeded 
2.5% of the network. This cut off is arbitrary but nearby val- 
ues do not appreciably change the resulting picture because 
above the phase transition the final size is bimodal: either 
spreading takes off and reaches a characteristic fraction of 
the network, or it fails. The network size is N = 10 5 and the 
resolution in r un( j is 0.01. See caption of Fig. [I] for further 
details. 



only type capable of receiving and transmitting an infec- 
tion. With these observations, and in setting (3 = 1 for 
simplicity, we obtain the coupled equations 



9 H = 

k± ,oo 



<»> fi 

fci ,oo \ 





ki,oc 





k\ ,oo 



(30) 



and 



Final Size of Infection 



We use the evolution equations given in Sec. IV 



_ to 

describe the growth of the spreading process on our 
example networks, with the main goal of determining the 
final size. We start with the equations for edge i nfee 
tion probabilities, , 



and 9 { ° ] 



k,t+V 



and we present their full model-specific 



Eqs. (13) and (14 1, 



forms in 



the 



Appendix. We let t - 
for fixed points, 0- 



> oo in these equations and solve 
and 9^ . Only two of the eight 



possible equations are coupled, those for and 9 



g(u) 



(Eqs. Al and A2), and thus they alone determine the 
final probabilities. As per our previous calculations, this 
collapse in equation number is because k\ nodes are the 





ki ,oo 



Tdir 



2(?l u) 

ki ,oo 



\ fei.oo/ 





k\ ,oo 





k\ ,oo 



(31) 

Solving both equations for 0v and equating the results 

leads to a cubic polynomial in 9? . One root is 9^ = 
and the others are solutions to 



= r dir (4 u) ) -r dir (r und + 2)4 u) 

\ A;i ,oo / k\ ,oo 

+ (r d ir + l)(r und + l) -2. (32) 



If 



We look for solutions for which < 9^ < 1. 

k\ ,oo 

(f d ir + 1) (T und + 1) < 2, we find the only feasible solu- 
tion is 9^ } = 0. When (r dir + l)(r und + 1) > 2 we find 

ki .oo 
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FIG. 3: Final size curves for example parameter choices for 
the toy model described in Sec. [V] with /3 = 1. Symbols 
indicate sample output from simulations with iV = 10 s and 
solid curves follow from Eqs.(l5|), |33l, |34|, |35|, and (36 1. 



For each value of r un( j and Tdi r in the plots A and B, we show 
a maximum of 10 values of (f)^ randomly chosen from 10 5 
individual simulations for which the final size exceeds cj>x > 
50 (or 0.05%). The overall fit between theory and simulation 
is excellent. See the captions of Fig. [I] and Fig. [2] for further 
simulation details. 



k\ ,oo 



= again but now also a non-trivial solution 

1 . • L (1 (T-dir + l)(7- U nd + l) 



(T un d+2) 



1 



. (r und + 2Y 



(33) 

As expected, the probability of an infected edge 8^ 

becomes nonzero as we move away from the phase tran- 
sition curve in (th t , T un d)-space, in agreement with our 
global spreading condition analysis, Eq. (22). Using our 

3 (u) 



expression for 



, we obtain expressions for the other 



non-zero edge infected probabilities: 

(1-T U 



(u) 

fel ,oo 



Tund(l-4 U) ) 
\ k\ ,oo / 



(34) 



k% ,oo 



= (1 - r dir ) 



9 W 

ki ,oo 



- 26i u) 

ki ,oo 



(35) 



-29^ 



\ ki,ooJ \ fci,oo k\,Q 



and 



9 (u) _ 

fe4 ,00 



= (1 - r und ) (8^ 

\ ki ,i 



fel ,00 



fel, oo fei,oo 



(36) 



Our last step is to use the above edge infection probabili- 
ties to compute the eventual fractional extent of a global 
spreading event 4>oo using Eq. (15) (with </>o — > 0). In 
Fig. [3j we compare output of our simulations with the 
model's version of Eq. ( 15 ), once again showing excellent 
agreement. 



VI. CONCLUDING REMARKS 

We have provided an extensive treatment of spread- 
ing on generalized random networks, accommodating 
a wide range of contagion processes from biological to 
social in nature. Our analysis is straightforward in that 
physical intuition is always at hand, and in no place 
have we resorted to more mathematical, less transpar- 
ent approaches, such as those employing generating func- 
tions. 

In closing, we note that if nodes are capable of recovery 
and reinfection, general calculations become considerably 
more difficult, particularly regarding the final extent of a 
spreading event, and this remains an open area of inves- 
tigation. 
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Appendix A: Final size calculations 

For the model described in Sec. [VJ we provide the spe- 
cific forms below for Eqs. (13) and (14) as used for the 
calculations in Sec. |V C| regarding the final size of spread- 
ing events. 
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